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ABSTRACT 

The formation timescale and final architecture of exoplanetary systems are closely related to the 
properties of the molecular disks from which they form. Observations of the spatial distribution and 
lifetime of the molecular gas at planet-forming radii (a < 10 AU) are important for understanding the 
formation and evolution of exoplanetary systems. Towards this end, we present the largest spectrally 
resolved survey of H2 emission around low-mass pre-main sequence stars compiled to date. We use a 
combination of new and archival far-ultraviolet spectra from the COS and STIS instruments on the 
Hubble Space Telescope to sample 34 T Tauri stars (27 actively accreting CTTSs and 7 non-accreting 
WTTSs) with ages ranging from ~ 1 - 10 Myr. We observe fluorescent H2 emission, excited by Lya 
photons, in 100% of the accreting sources, including all of the transitional disks in our sample (CS 
Cha, DM Tau, GM Aur, UX Tau A, LkCal5, HD 135344B and TW Hya). The spatial distribution 
of the emitting gas is inferred from spectrally resolved H2 line profiles. Some of the emitting gas is 
produced in outflowing material, but the majority of H2 emission appears to originate in a rotating 
disk. For the disk-dominated targets, the H2 emission originates predominately at a < 3 AU. The 
emission line-widths and inner molecular radii are found to be roughly consistent with those measured 
from mid-IR CO spectra. 

Subject headings: protoplanetary disks — stars: pre-main sequence — ultraviolet: planetary systems 



1. INTRODUCTION 

The lifetime, spatial distribution, and composition of 
gas and dust in the inner ~ 10 AU of young (age < 10 
Myr) circumstellar disks are important components for 
understanding of the formation and evolution of extra- 
solar planetary systems. The formation of giant planet 
cores and their accretion of gaseous envelopes occurs on 
timescales similar to the lifetimes of the disks around 
Classical T Tauri Stars (CTTSs; 10 6 - 10 7 yr). The 
cores of giant planets are th ought to be compri sed of 
coagulations of dust grains (jHavashi et al.1 IT985). and 
the majority of observational work on the lifetime of in- 
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ner disks has come from photometric and spectroscopic 
studies of their dus t (jHaisch et al.ll200fllHernandez et al.1 
120071: iWvattl 120081 ). Dust in protoplanetary disks is ob- 
served as mid- and far-IR excess flux produced by warm 
grains (e.g., Furlan et al. 2006, 2009; Evans et al. 2009; 
Luhman et al. 2010). The IR spectral energy distribu- 
tions (SEDs) of protoplanetary disks are sensitive to the 
radial distribution of dust in the disk; this dependence 
has led to the discovery of a class of "transitional" sys- 
tems, whose SEDs indicate that gaps of a few tenths to 
tens of AU have been opened in their inner disks (Strom 
et al. 1989; Calvet et al. 2002, 2005; Espaillat et al. 
2007; and see the review by Williams & Cieza 2011). 
The physical process by which the inner disk is cleared 
is not yet establishe d. Possible mechanisms includ- 
ing p hotoevaporation (|Alexander et a l. 2006; Gort i et al.l 
2009) and dynamical clearing by exoplanetary sys - 
tems (jRice et al.ll2003HDodson-Robinson k Salvldl20ill ). 
pos sibly aided by the magn etorotational instabil- 
ity (|Chiang fc Murrav-C lav 2007), can reproduce certain 
transitional disk observations. 

The lifetime and spatial extent of the ga s disk de- 
term ine the final mass of giant planets (Ida fc Lin! 
120041 ) and the final architecture of an exoplanetary 
syste m, as disk g a s regulates type-II planetary migra- 
tion flWardJ 119971 : lArmitage etail l2002t iTrilling et all 
2002). Because the migration timescale is sensitive 
to the specifics of the disk s urface density distribu- 
tion and dissipation timescale (jArmitage! 120071 ). obser- 
vations of gas-rich systems with age < 10 Myr can 
provide important constraints on models of the evo- 
lution of exoplanetary systems. Significant observa- 
tional effort has been devoted to the study of inner 
disk gas in recent years, including ground-based mid- 
IR spectroscopy of CO and [Ne II] (jNaiita et all 120031: 
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Pascucci etaLl 120071: IHerczeg et"al) 120071: iNaiita et all 
20091 iBast et al.l I201lt iSacco et al.N2012D. spectroastro- 



metr ic observati ons of CO (iPo ntoppid an et al.l [2008a, 
120111 ) and [Ne II] (jPascucci et al. 1120 111), and Spitzer-IRS 
observations of H?Q and organics dCarr k Naiital 120081 : 
iSalvk et al.ll200ll2011bl: iCarr k Najitall201lD . There is 
growing evidence that remnant gas disk survival is com- 
mon inside the dust hole in transitional disks (e.g., Salyk 
et al. 2009, 2011), suggesting that planetary migration 
may continue after the dust disk has dispersed. In tran- 
sitional systems with minimal inner disk dust, observa- 
tions of active accretion also provide i ndirect evidence for 
the presence of a remnant gas disk (INaiita et al.l l2007t 
IKim et al.lf2q09t IMerin et al.ll2O10HFedele et al.ll2O10fL 

Many previous studies of inner disk gas have used trace 
species to infer the presence of molecular hydrogen (H2), 
the primary constituent of protoplanetary disks and gas 
giant planets. The homonuclear nature of H2 means that 
rovibrational transitions are dipole forbidden, with weak 
quadrupole transitions that have large energy spacings. 
This makes direct detection of H2 challenging at near- 
and mid-IR wavelengths (Pascucci et al. 2006; Carmona 
et al. 2008; but see also Bary et al. 2008). However, 
H9 can be observed in the far-UV (912 - 1650 A) band- 
pass, where the strong dipole-allowed electronic transi- 
tion spectrum is prim arily photo-excited ( "pumped" ) by 
stella r Lya photons (|Ardila et al.l 120021 : IHerczeg et al.l 
12001 . The Lya -pumping route proceeds primarily by 
absorption out o f the second excited vibrational level 
(v = 2) of H 2 (|Shulll 119781 ). which implies that the 
molecules reside in a hot (T(H 2 ) > 2000 K) disk sur face 
at semi- major axes a < 10 AU (IHerczeg et al. 2004), or 
in ex tended outflows (| Walter et al.l 120031 : iSaucedo et al.1 
2003). An analysis of the spectral line profiles can dis- 
tinguish between these origins, enabling one to identify 
and characterize emission from the molecular disk given 
sufficient spectral resolution and spectroscopic sensitiv- 
ity. Assuming that disk molecules are in Keplerian orbit 
around their central star, line broadening due to orbital 
motion dominates the profile in moderate-to-high incli- 
nation systems. H2 velocity widths can therefore be used 
to infer the spatial distribution of the gas on the disk sur- 
face. 

In this work, we present the most sensitive survey of 
spectrally resolved H2 emission in protoplanetary disks 
obtained to date. Previous spectral sur veys of TTSs 
with the International Ultravi olet Explorer (|Valenti et al.1 
2000; Johns-Krull eTaLl l2000l) and th e various ultraviolet 



spect r ographs on HST (lArdila et al.ll2002l: IHerczeg et al.l 
[20061: llnglebv et ail 120091: lYang et al.l I2012D have been 
carried out at either lower spectroscopic sensitivity or 
resolution. In this study, we take advantage of the high 
sensitivity, low instrumental background, and moderate 
spectral resolution of the Hubble Space Telescope-Cosmic 
Origins Spectrograph to greatly expand the number of 
targets available for detailed UV studies. In §2, we de- 
scribe the targets and the HST observations. The anal- 
ysis performed to characterize the H 2 luminosities and 
spatial distributions are described in §3. We present in 
§4 a discussion of H 2 line profiles, considering outflow and 
disk origins for the emitting gas, and an estimate of the 
fraction of stellar Lya re-processed by circumstellar H2. 
§4 also presents the time-evolution of the amount and 
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Fig. 1.— Examples of COS spectra (1170 - 1690 A) for a range of 
gas and dust disk parameters. From top to bottom: The primordial 
disk target AA Tau [offset by +150 (x 10~ 15 erg cm -2 s _1 A -1 )], 
the pre-transitional disk V4046 Sgr [flux divided by 4 and offset by 
+75 (x 10 -15 erg cm -2 s _1 A -1 )], the transitional disk system 
GM Aur, and the gas-poor WTTS LkCal9 [offset by -50 (x 10~ 15 
erg cm -2 s — 1 A -1 )]. The spectra have been binned by one spectral 
resolution element (7 pixels) for display. Except for the atomic lines 
identified in the spectrum of LkCal9, most of the emission lines 
in the spectra of the other three stars are fluorescent H2 emission 
lines pumped by Lya. 



location of the H2 gas, suggesting that the ^-emitting 
gas both dissipates and moves towards larger orbital radii 
over the interval from 10 6 - 10 7 yr. A brief summary of 
the results from the molecular survey are presented in 
§5- 

2. TARGET SAMPLE AND OBSERVATIONS 

2.1. Target Sample 

The goal of this observational survey is to span a range 
of ages, mass accretion rates, and star-forming environ- 
ments in order to better understand the global properties 
of H2 emission in protoplanetary environments. The tar- 
gets mainly belong to the Taurus- Auriga, r\ Chamaeleon- 
tis, TW Hya, and Chamaeleon I star-forming regions, as 
well as individual targets in other associations and iso- 
lated systems. Potential sources of uncertainty in ana- 
lyzing a diverse population of targets are systematic ef- 
fects based on different methods used to derive system 
parameters in the literature. In order to mitigate the 
effects of systematics on this study of molecular disks, 
our approach was to adopt system parameters from pa- 
pers where similar techniques were used to derive prop- 
erties such as ages, extinctions, stellar masses, and in- 
clinations. The adopted target parameters are given 
in Table 1. Where possible, 1) inclinations were taken 
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TABLE 1 
HST Target List 



Target a 


- 1 V 


logi n ( Age) 


M* 


M 


/. 


HST a 


Ref. c 






(y rs ) 




{lv ivi0 yr ) 


y°) 


PID 




AA Tau 


0.50 


6.38 ± 0.20 


0.80 


0.33 


75 


11616 


2,15,28 


AK Sco 


0.5 


7.24 ± 0.24 


1.35 


0.09 


68 


11616-S 


3,24 


BP Tau 


0.50 


5.94 d= 0.29 


0.73 


2.88 


30 


12036 


1,15,30 


CS Cha 


0.8 


6.39 ± 0.09 


1.05 


1.20 


60 


11616 


4,25 


CV Cha 


1.67 


6.70 ± 0.10 


2.00 


3.16 


35 


11616-s 


5,26 


DE Tau 


0.60 


5.82 ± 0.20 


0.59 


2.64 


35 


11616 


1,15,29 


DF Tau A 


0.60 


6.27 ± 0.53 


0.19 


17.7 


85 


11533 


15,29 


DK Tau A 


0.80 


6.17 ± 0.22 


0.71 


3.79 


50 


11616 


1,15,29 


DM Tau 


0.0 


6.56 ± 0.20 


0.50 


0.29 


35 


11616 


2,16,21 


DN Tau 


1.90 


6.04 ± 0.20 


0.60 


0.35 


28 


11616 


2,15,21,31 


DR Tau 


3.20 


6.18 ± 0.20 


0.80 


3.16 


72 


11616 


2,17,28 


GM Aur 


0.10 


6.86 ± 0.20 


1.20 


0.96 


55 


11616 


2,15,21 


HD 104237 


0.70 


6.30 ± 0.30 


2.50 


3.50 


18 


11616-S 


7,19,40 


HD 135344B 


0.30 


6.90 ± 0.30 


1.60 


0.54 


14 


11828 


8,19,36 


HN Tau A 


0.5 


6.27 ± 0.27 


0.85 


0.13 


>40 


11616 


1,15,22 


IP Tau 


0.20 


6.37 ± 0.24 


0.68 


0.08 


60 


11616 


1,15,32 


LkCa 15 


0.60 


6.35 ± 0.26 


0.85 


0.13 


49 


11616 


1,16,21 


RECX 11 


0.0 


6.60 ± 0.20 


0.80 


0.03 


70 


11616 


9,20,45 


RECX 15 


0.0 


6.78 ± 0.08 


0.40 


0.10 


60 


11616 


10,20 


RU Lupi 


0.07 


6.39 ± 0.09 


0.80 


3.00 


24 


12036 


11,18,34 


RW Aur A 


1.6 


5.85 ± 0.53 


1.40 


3.16 


77 


11616 


1,23,35 


SU Aur 


0.9 


6.39 ± 0.21 


2.30 


0.45 


62 


11616 


1,17,37 


SZ 102 


1.13 


6.15 ± 0.15 


0.75 


0.08 


10 


11616 


12,27,38 


TW Hya 


0.0 


7.00 ± 0.40 


0.60 


0.02 


4 


8041-5 


13,18,36 


UX Tau A 


0.20 


6.10 ± 0.30 


1.30 


1.00 


35 


11616 


1,21 


V4046 Sgr 


0.0 


6.90 ± 0.12 


0.86+0.69 


1.30 


36 


11533 


14,22,39 


V836 Tau 


1.70 


6.26 ± 0.26 


0.75 


0.01 


65 


11616 


1,18,41 


HBO 427 


0.00 


6.64 ± 0.14 


0.7 






11616 


6 


LkCa 19 


0.00 


6.84 ± 0.38 


1.35 






11616 


1 


LkCa 4 


0.69 


6.43 ± 0.25 


0.77 






11616 


6 


RECX 1 


0.00 


6.78 ± 0.00 


0.90 






11616 


12 


TWA 13A 


0.00 


6.90 ± 0.12 


0.32 






12361 


43 


TWA 13B 


0.00 


6.90 ± 0.12 


0.38 






12361 


13 


TWA 7 


0.00 


6.39 ± 0.39 


0.55 






11616 


44 



a Targets in the upper group are CTTSs, targets in the lower group are WTTSs. 
b Program IDs marked — S indicate that STIS observations were used 

c (1) IKraus fc Hillenbrand! ( 120091') ; (2) IRicci et all IpOlOl'). aee uncertain ties are assumed to be ± .20; (3) lAlencar et alJ [|20"03T) ; 
(4) Lawson et al. (1996); (5) Sicss ct al. (2000); (6) Bcrtout et al. (2007); (7) Feigclson et al. (2003a); (8) van Boekel et al. (2005); 
(9) Lawson ct al. (2001); (10) Ramsav Howat & Greaves (2007); (11) Hcrczcg et al/ ( 2005a); (12) Comcron & Fernandez ( 2010); (13) 
[ Webb et al.l H19991); (14 1 Quast ct al. (2000a); (15) Gullbring ct al. (1998); (16) Hartman n et al.l fl998i1: (17) IGuIIbring et al.) 112000); (18) 
Herczeg & Hillenbrand (2008); (19) Garcia Lopez ct al. (2006); (20) Lawson ct al. (2004); (21) Andrews ct al. (2011); (22) France et al. 
(2011b); (23) White & Ghcz (2001); (24) Gomez dc Castro (2009); (25) Espaillat et al. (2007a); (26) Hussain ct al. (2009); (27) 
Comcron ct al. (2003a); (28) Andrews & Williams (2007); (29) Johns-Krull & Valcnti (2001); (30) Simon ct al. (2000); (31) Muzcrollc ct al. 
(12003 1 ) ; (32T[Espaillat ct al. (2010T); ( 34) IStempels et al.1 J2007f ); (3 5)IEisner et al.1 (120071) ; (36) Pontoppidan ct al. (2008b); (37)IAkeson eTaL 
(2002) : (38)ICoffev et a l. (2004); (39) Rodri guez et al.l 1120101 ) : ( 40)IGradv et aljlfeoolh (41) INaiita et al.l [|2008r )~ (42) IInglebv et al.l H2011bl ): 
(43) IPlavchan et al.l (120091) ; (44) INeuhauser et al.l (12000) 7745) llnglebv et al.l 1 1201 lal ) 



from sub-mm/IR interferometric studies, 2) ages, stel- 
lar masses, and extinctions were derived from pre-main- 
sequence stellar evolutionary tracks, and 3) mass accre- 
tion rates were derived from measurements of the accre- 
tion luminosity. Of our 34 targets, 27 are considered 
CTTSs while 7 do not show evidence for active accre- 
tion or a gas-rich circumstellar disk and are classified as 
Weak-lined T Tauri Star (WTTSs). These populations 
are separated as the upper and lower groups, respectively, 
of target stars listed in Table 1. 7 of the CTTSs in our 
sample are considered transitional systems (CS Cha, DM 
Tau, GM Aur, UX Tau A, LkCal5, HD 135344B and 
TW Hya), with mid-IR SEDs indicating that a gap has 
opened in their inner dust disks. This sample includes 
several of the best-studied transitional disks in the lit- 
erature. References are listed in Table 1, and a more 
detailed des cription of a subsample o f these targets can 
be found in iSchindhelm et ail (|2012a[) . 
The assumed distances are not critical to the results 



presented here, but they do impact the comparison of 
the H2, Lya, and C IV luminosities presented in §4. For 
the Taurus- Auriga targets (AA Tau, BP Tau, DE Tau, 
DF Tau, DK Tau, DM Tau, DN Tau, DR Tau, GM Aur, 
HN Tau, IP Tau, LkCal5, RW Aur, SU Aur, UX Tau A, 
V836 Tau, HBC 427 (V397 Tau), LkCal9, LkCa4), we 
assumed d — 140 pc (Elias 1978; Kenyon & Hartmann 
1995; and see also the VLBA work presented by Loinard 
et al. 2007); for the 77 Cha targets (RECX-11, RECX-15 , 
RECX-1), we assumed d = 97 pc (|Mamaiek et al.lll999l) : 
for the TW Hya association targets (TW H ya, TWA13A 
TWA13B, TWA7), we assumed d = 54 pc ()van Leeuwenl 
[20071) . and for Cha maeleon I (CS Cha, CV Cha), we 
assumed d — 160 pc (|Luhmanl [20041 . Other objects are 
V4046 Sgr (d = 83 pc; Quast et al. 2000), SZ 102 (in 
Lupus 3, d = 200 pc; Comeron et al. 2003), HD104237 (a 
member of the e Cha group, d = 116 pc; Feigelson et al. 
2003), AK Sco (d = 103 pc; van Leeuwen 2007), and RU 
Lupi (d — 121 pc; van Leeuwen et al. 2007). Six of the 
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TABLE 2 
Selected H2 Emission Lines. 



1440 1450 o 1460 
Wavelength (A) 



1470 



Fig. 2.— The 1430 - 1470 A spectral region for the gas-rich 
targets plotted in Figure 1. All of the strong spectral features 
in this bandpass are emission lines from Lya-pumped fluorescent 
H2. Emission lines used in this analysis are marked with blue di- 
amonds and several bright features are labeled. Objects are plot- 
ted in order of decreasing ncar-IR dust excesses: AA Tau (pri- 
mordial, ni 3 _3i = -0.51; Furlan et al. 2009), V4046 Sgr (pre- 
trans itional) has a sub-AU sca le hole in the inner disk dust distribu- 
tion ( Jensen & Mathicu 1997), and GM Aur (transitional, ni3— 31 
= 1.7 5; Furlan et al. 20 09) has a ~ 24 AU hole in the inner dust 
disk {Calvct et al. 200||). The spectra are binned to » 1/2 of a 
spectral resolution element (3 pixels), and representative error bars 
are shown overplotted. The H2 spectra are qualitatively similar, 
independent of the inner disk dust properties. 



nineteen sources in Taurus- Auriga are known multi-star 
systems (see, e.g., Kraus et al. 2012) and V4046 Sgr is 
known to be a short-period binary system (Qu ast et al.l 
l2000bD . 

2.2. Observations 

Our sample of 34 T Tauri stars was assembled from 
new and archival observations with HST-COS and - 
STIS. The majority of the targets were observed as part 
of the DAO of Tau guest observing program (PID 11616; 
PI - G. Herczeg) and the COS Guaranteed Time Observ- 
ing program (PIDs 11533 and 12036; PI - J. Green). Ad- 
ditional observations of the transitional disk HD135344B 
and weak-lined systems TWA13A and TWA13B (PIDs 
11828 and 12361; PI - A. Brown) are presented. A sub- 
set of the H2 survey obser vations have been p re sented 
previously in the lite rature dFrance et al.ll2011bllaL l2012t 
Indebv et afll2011bt lYang et al.H201lHSchindhelm et al.l 



2012a|). Finally, we have included archival STIS observa- 



tions of the well-studied CTTS T W Hva jlHerczeg et al.l 
2002), obtained through StarCAT (|Avresll2010| ) 

Most of the targets were observed with the medium- 
resolution far-UV modes of COS (G130M and G160M; 
Green et al. 2012). These observations were acquired 
between 2009 December and 2011 September. Multiple 
central wavelength settings at several focal-plane split 
positions were used to create continuous far-UV spectra 
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a Transitions are for the B 1 ^^ - H2 band system. 

The branching ratio is the ratio of the line transition probabil- 
ity to the total transition probability out of state [v ,J ], Bmn — 

A n ','j'-^".j" 
— ^JZrj, — 

c Blended with Si iv Xi ab = 1402.77 A in some targets. 

from w 1150 - 1750 A and mitigate the effects of fixed 
pattern noise. These modes provide a point-source reso- 
lution of Af w 17 km s"" 1 wit h 7 pixels per resolution el- 
ement (jOsterman et al.ll201lD . The data were smoothed 
by three pixels for analysis. The total far-UV exposure 
times were between two and four orbits per target, de- 
pending on the intrinsic luminosity and the interstellar 
plus circumstellar reddening on the sightline. The one- 
dimensional spectra produced by the COS calibration 
pipeline, CALCOS, were aligned and co added using the 
custom software procedure described by iDanforth et al.l 
(2010). The full far-UV spectra of three CTTSs (AA 
Tau, V4046 Sgr, and GM Aur) and one WTTS (LkCal9) 
are displayed in Figure 1, and a 40 A blow-up of the 
CTTSs is shown in Figure 2. 

Targets that exceeded the COS bright-object limit 
were observed with STIS in the medium-resolution 
echelle mode. We used the E140M mode (Av 8 km 
s^ 1 ) mode through the 0.2" x 0.2" slit for exposure times 
between two and three orbits per object. The far-UV 
STIS spectra were combined using the STIS echelle soft- 
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Fig. 3a.— Velocity profiles for the H 2 S 1 Sj - X 1 ^ (1 - 7) R(3) (A !ab = 1489.57 A) emission line in all targets, plotted over the ± 
200 km s _1 interval. The data plotted here are smoothed to 3 pixels (fs 0.5 spectral resolution elements), with representative error bars 
shown in orange. The green dash-dotted line is plotted at the stellar radial velocity, when known. The target labels are the target names 
identified in the headers of the HST spectroscopic observations. 



ware developed for the StarCAT catalog (T. Ayres - pri- 
vate communication; Ayres 2010). 

3. ANALYSIS 

We observe fluorescent H2 emission from all of the 27 
CTTSs in our sample, and no H2 emission from the 7 
WTTS targets. The number of observed fluorescent pro- 
gressions varies significantly across the sample, and we 
present our measurements of the total H 2 fluxes below. 
We detect strong fluorescent emission in all of the transi- 
tional objects in our sample (CS Cha, DM Tau, GM Aur, 
UX Tau A, LkCal5, HD 135344B and TW Hya; §4.3). 

The fluorescent H2 lines observed in the CTTS sample 
can be used to determine the relative amount of H 2 in 
the circumstellar environment and to constrain its spa- 
tial distribution. For our line-profile analysis, we fo- 
cus on the measurement of two progressions {[v ,J ] = 
[1,7] and [1,4]) detected in all of the CTTS targetQ 

12 The quantum numbers v and J denote the vibrational and ro- 
tational quantum numbers in the ground electronic state (X 1 S^") ) 
the numbers v and J characterize the H2 in the excited (B 1 Sj) 
electronic state, and the numbers v and J are the rovibrational 



These emission lines are pumped through the (1 - 2)R(6) 
A ;ab 1215.73 A and (1 - 2)P(5) A Zab 1216.07 A transitions, 
respectively. The absorbing transitions are within +15 - 
+100 km s _1 of Lya line-center. The signal-to-noise ra- 
tios (S/N) per resolution element are typically between 
5 and 40 in the brightest fluorescent H2 emission lines 
for our CTTS targets. The (1 - 7)R(3) A ;ab = 1489.57 A 
transition is relatively free from spectral contamination 
and is displayed for all targets in Figure 3 (a - c) . When 
available, the stellar radial velocities are indicated in Fig- 
ure 3 with green dash-dotted lines. Within the ~ 15 km 
s _1 wavelength solution accuracy of COS, most of the H2 
progressions are consistent with the stellar velocity. In 
several cases the H2 lines appear to have line wings ex- 
tending to negative velocities or statistically significant 
differences in the velocity widths of different H2 progres- 
sions. RW Aur is an e xtreme example of th is behavior. 
lArdila et all (|2002f ) and lHerczeg et al.l (|2006D have noted 
the presence of blue-shifted H2 emission in their sample 

levels in the electronic ground state following the fluorescent emis- 
sion. Absorption lines are described by (v — v) and emission lines 
by (v - V ). 
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of T Tauri stars observed with GHRS and STIS. The 
presence of blue-shifted emission creates additional un- 
certainty in the measured line-widths at the resolution 
of our spectra, and we present a discussion of outflow 
signatures observed in our sample in §4.2.1. 

The high S/N of the COS data means that we can re- 
strict the analysis to the brightest lines from the progres- 
sions of interest. Systematics were minimized by focusing 
on emission lines with wavelengths 1395 < A < 1640 A. 
This choice is optimal because 1) fluorescent transitions 
cascading to vibrational levels v > 5 do not suffer sig- 
nificant self- absorption before escaping the circumstel- 
lar environment (e.g., Figure 7 of Herczeg et al. 2004), 
enabling more robust flux measurements, 2) the non- 
Gaussian wi ngs of the COS instrumental line-spread 
function (LSF[3) contain a smaller fraction of the power 
at A > 1400 A, 3) the strongest lines from the [1,4] and 
[1,7] progressions are in this range, and 4) the relative 
correction for interstellar reddening between the differ- 
ent fluorescent emission lines is minimized. The first and 

13 The LSF experiences a wavelength dependent 
non-Gaussianity due to the introduction of mid- 
frequency wave-front errors produced by the polish- 
ing errors on the HST primary and se condary mirrors; 
http : //www . stsci . edu/hst/cos/documents/isrs/ 



latter two arguments apply equally to the STIS data. 



3.1. Line Fluxes and the H2 Progression Luminosity 

The H2 emission lines were fit with an interactive 
multi-Gaussian IDL line-fitting code optimized for COS 
emission line spectra. This code assumes a Gaussian line- 
shape convolved with the wavelength dependent LSF, 
then uses the MPFIT routine to minimize \ 2 between 
the fit and data ([Markwardt 2009|). A second order poly- 
nomial background, the Gaussian amplitudes, and the 
Gaussian FWHMs for each component are free parame- 
ters. The parameters of the underlying Gaussian emis- 
sion lines are returned to the user. The smaller STIS 
aperture does not sample the broad wings of the HST 
LSF, therefore unconvolved Gaussians were used for the 
targets observed with STIS. We chose 12 progressions 
with absorbing transitions that spanned the width of the 
observed Lya profiles in most targets, 1213.3 < A a i, s (H 2 ) 
< 1219.1 A (jSchindhelm et al.ll2012bD . We fit the bright- 
est unblended lines from 12 progressions in the 1395 - 
1640 A bandpass, and used these values to determine 
the total flux from each progression. 

A list of the selected lines is given in Table 2. The total 
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flux from a given progression is given by 

'-(*>- *E(lb) <" 

where F mn is the reddening corrected, integrated H 2 
emission line flux from rovibrational state m (= [v , J ]) 
in the -B 1 ^ electronic state to n (= [v ,J ]) in the 
ground electronic state, X 1 E+ B mn is the branching 
ratio between these two states, and N is the number of 
emission lines measured from a given progression. The 
measurement errors are typically small, so we take the 
flux error to be the standard deviation of the individ- 
ual measurements of F m (H.2). The dominant systematic 
error on the measured H2 flux is the correction for inter- 
stellar reddening; we do not attempt to account for this 
uncertainty in the flux and luminosity errors presented 
below. See the ApJ version of this paper (or contact the 
authors) for reddening-corrected progression fluxes for all 
of the gas-rich targets. The total progression luminosity 
is then L m (H2) = (47r<i 2 )F TO (H2). Upper limits on the 
H2 emission line fluxes in the gas-depleted targets were 
determined from the standard deviation in a ± 50 km 
s _1 region surrounding the laboratory wavelength of the 
transition. The total fluorescent H2 luminosity is then 



taken as the sum of all 12 of the fluorescent progressions. 
This prescription should account for > 80 % of the to- 
tal Lya-pumped H2 emission from our targets, although 
the exact fraction of the total measured H2 flux will de- 
pend on the local Lya line-profile (Herczeg et al. 2004; 
Schindhelm et al. 2012b). 

Depending on the geometry and spatial distribution of 
the absorbing molecular layer, the H2 absorption lines 
may be optically thin (or have optical depths of a few) 
in some targets. In this case, the emitted H2 luminos- 
ity will be directly proportio nal to the number of Lya 
pumping photons received. iSchindhelm et all {2012b) 
have presented total incident Lya fluxes for 14 of the 
targets in our sample. Using a line-profile reconstruc- 
tion technique that takes into account the 12 H2 pro- 
gressions described above, they simultaneously fit the 
neutral hydrogen outflow, the H2 column density, and 
the H2 temperature (previous examples of H2-based Lya 
reconstructions are described by Wood et al. 2002; Wood 
k Karvoska 2004; Herczeg et al. 2004). The amount of 
the stellar Lya that is redistributed by H2 can be de- 
termined by dividing the total H 2 flux, i^(H 2 ), by the 
total Lya flux, F(~H.2)/F(Lya). The H2 luminosity can 
also be compared to accretion indicators in our data set, 
such as C IV. Emission from the C IV AA1548,1550 A 
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Fig. 4.— Expanded view of the H 2 (1-6)P(5) (\ lab = 1446.12A) 
emission line, showing the velocity width beyond the COS line- 
spread function (shown as the gray dash-dotted line) for the three 
example targets shown in Figure 2. The A A Tau and V4046 Sgr 
spectra have been normalized to the peak flux of the GM Aur H2 
emission and shifted to the centroid velocity of the H2 emission. 
The emission lines are spectrally resolved by COS. 

resonance doublet in excess of a baseline magnetospheric 
level has been shown to correlate with the mass accre- 
tion rate (IJohns-Krull et al.l2000HYang et al.l20ll) . The 
C IV flux can be measured directly in the same data set 
from which we measure FmQUz), mitigating complica- 
tions associated with short- and long-baseline time vari- 
ability. Due to the non-Gaussian appearance of many of 
the C IV emission profiles, we measured F(C IV) by in- 
tegrating the reddening corrected spectra over (1547.5 - 
1553.5 A) and subtracted the continuum from an adja- 
cent, line- free portion of the spectrum. 

3.2. H2 Line Widths and the Average H2 Radius, (RH2) 

Figure 4 shows the H 2 (1 - 6)P(5) (A; ?6 = 1446.12 A) 
profiles of the three example spectra displayed in Fig- 
ure 2, with the 1450 A HST+COS LSF overplotted as 
the gray dash-dot line. These spectra are typical of the 
CTTS sample and one observes that the emission lines 
are spectrally resolved. We display the basic H2 and C IV 
line parameter observations for the [1,7] and [1,4] pro- 
gressions in Figure 5. The H2 emission line FWHMs are 
uncorrelated with the stellar mass, i^ m (H 2 ), or F(C IV). 

Kinematic broadening dominates the observed H2 line 
profiles (see §4.2). The thermal broadening of the emis- 
sion lines is approximately 4.5 km s _1 at the nominal 
2500 K H2 layer, significant additional broadening would 
require temperatures in exc ess of the ~ 4500 K disso- 
ciation temperature of H 2 (jLepp fc Shulll [l983() . If we 
further assume that any turbulence in the disks is sub- 
sonic, then the maximum turbulent velocity will be no 
larger than a few km s _1 . Therefore, velocity broaden- 
ing due to bulk motions and Keplerian rotation domi- 
nate the observed line shapes when the FWHM of the 
emission line is greater than the 17 km s spectral reso- 
lution of COS. For the case of H2 in a circumstellar disk, 
we define a simple metric to characterize the average H2 



radius, (R112), 

where M* is the stellar mass, i is the inclination an- 
gle, and FWHM m is the mean of the Gaussian FWHMs 
for a given progression m. This definition of the aver- 
age mol ecular radius is an a logous to the UV-CO radius 
used by iSchindhelm et~aLl (|2012af ). When possible, we 
use disk inclinations derived from the sub-mm dust con- 
tinuum obs ervations presented by lAndrews fc Williams! 
(|2007| ) and lAndrews et al.1 (|2011l ). In principle, the er- 
ror on the average H2 radius should include uncertainties 
on the stellar mass and disk inclinations, however these 
uncertainties are not available for all targets. Therefore, 
the quoted error on {Rh2) only includes measurement 
uncertainties from the H2 line fitting. 

In Table 3 we present {Rh2) for the [v , J ] = [1,7] 
progression (= (Rh2)[i,7])- We focus on the [1,7] pro- 
gression due to its proximity to the Lya line-center. This 
progression should be one of the most readily observable 
in weakly accreting systems because as the accretion- 
powered Lya flux decreases, the flux in the wings of 
Lya may be insufficient to excite a detectable level of 
H2 emission. However, the [1,7] progression should con- 
tinue to be observable even in systems with narrow Lya 
profiles (e.g., France et al. 2010). While we concentrate 
on the (i?^f2)[i,7] , the (Rh2)[ia] distribution is qualita- 
tively similar. We do not compute (Rh2) for targets with 
inclination angles < 15°; the small radial component of 
the H2 velocity vector makes the derived radii very sen- 
sitive to uncertainties in the disk geometry. We also do 
not compute (Rh2) for the only unambiguously outflow- 
dominated source in the sample, RW Aur (§4.2.1). The 
impact of a second, weaker outflow emission component 
will bias our results towards smaller H2 radii in some tar- 
gets, but a single component dominates the majority of 
our sources, therefore we adopt a single emission compo- 
nent in order to facilitate a uniform spectral analysis. 

A knowledge of the inner disk radius, Rim is impor- 
tant for understanding the star-disk interaction. In or- 
der to avoid complications from outflows or blending 
from adjacent weaker H2 transitions, we do not mea- 
sure i?m(H2) from half- width zero-intensities (HWZI, 
e.g., Brittain et al. 2009 ). Instead, we adop t the defini- 
tion of Rin suggested bv lSalvk et al.l (|2011a[ ) for mid-IR 
CO emission from the inner disk, where is the Ke- 
plerian semi-major axis corresponding to 1.7 x HWHM 
of the Gaussian line fit. This choice leads to the rela- 
tion between the average H2 radius and the inner H2 
radius, (i?H2)[i,7] = 2 -89 i? m (H 2 ). i?i„(H 2 ) is given with 
(Rh2}[i,7] in Table 3. 

4. DISCUSSION 

There are numerous indicators of the gas and dust con- 
tent of a young protoplanetary system. Three important 
observables are the warm dust content of the inner disk, 
the presence of circumstellar gas, and signs of active ac- 
cretion. Our HST observations provide measurements 
of the last two, while the first has been extensively stud- 
ied in the IR. The combination of low spectral resolution 
and large instrumental backgrounds that have compli- 
cated the detection and analysis of H 2 in previous UV 
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Fig. 5. — The average H2 FWHM is compared with the stellar mass, the progression flux (_F m (H2)), and the C IV flux (F(C IV)) 
for two progressions: [1,7] (pumped through the (1 - 2) R(6) \i a b = 1215. 73A transition) and [1,4] (pumped through the (1 — 2) P(5) 
^lab = 1216. 07A transition). The H2 emission line FWHMs are uncorrelated with the stellar mass, i*Vn(H2), or F(C IV). 



surveys has largely been remedied with the installation 
of COS. The large transition probabilities of the H2 elec- 
tronic band systems and the lack of photospheric emis- 
sion at A < 1700 A in low-mass stars make fluorescent 
H2 one of the most sensitive indicators for the presence 
of molecular gas in the inner ~ 10 AU of young circum- 
stellar disks. 

Accretion shocks are a significant source of hot 
gas in accreting systems, observed as UV and X-ray 
emission lines in excess of what can be attributed 
to magnetospheric activity alone dCalvet fc Gullbrind 
[l998t I Johns-Krull etaTII^OOOt IGiinther fc Schmittl l2008^ 
Specifically, excess emission from neutral hydrogen (line 
formation temperature Tf orm ~ 10 4 K; observed as Lya 
and Ha emission) and the C 3+ ion (Tf orm ~ 10 5 K; ob- 
served through the AA 1548, 1550 A C IV resonance dou- 
blet) correlate well with both the m ass accretion rate 
and t he H2 emission from the system (jJohns-Krull et ahl 
2000). This supports a symbiotic picture where gas- rich 
disks provide fuel for active accretion, and that accretion 
dominates the production of the Lya photons that make 
the H2 disk detectable. 

The total H2 luminosity is compared with the Lya and 
C IV luminosities in Figure 6. The general trend follows 



the expected relation that systems with larger L(Lya) 
and L(C IV) are actively accreting gas-rich disks. Lya 
fluxes are only available for about half of our sample (see 
Schindhelm et al. 2012b), and we interpolate (or extrap- 
olate) the strong correlation between F(Lya) and F(H.2) 
to estimate L(Lya) for the remaining objects. The direct 
measurements are shown as black squares in Figure 6, the 
interpolated (extrapolated) values are shown in green, 
and upper limits are indicated in red. The correlation 
between L(H2) and L(C IV) has a spread of ~ 1 - 1.5 
orders of magnitude in £(H2) at a given C IV luminos- 
ity. This spread may be partially due to uncertainties 
in the distance and reddening correction used to derive 
the luminosities, disks/outflows with differing amounts of 
molecular gas, and intrinsic variations in the C IV flux. 
In the following subsections, we combine our molecular 
and atomic tracers to constrain the properties of H2 in 
the circumstellar environments of these systems. 

4.1. Evolution of the Hi Luminosity and Reprocessing 
of the Lya Radiation Field 

The dissipation of inner dust disks (at a < 1AU) is 
thought to be mostly complete by w 6 Myr (e.g., Haisch 
et al. 2001; Wyatt 2008 and references therein). If the 
molecular gas and dust are coupled, we would expect 
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Fig. 6. — A comparison of the hot gas [logio(Ty orm ) ~ 4 — 5] and molecular gas emission from the TTSs studied here. The Lya 
and C IV emission is produced in the protostellar atmosphere through a combination of magnetic activity and magnetically funneled 
accretion. The H2 resides in the circumstcllar disk and in some cases, extended outflows, excited by Lya photons. The H2 luminosity is 
the summed luminosity of the 12 progressions listed in Table 2, representing > 80% of the total H2 luminosity (depending on the specific 
Lya pumping profile). Targets with H2 upper limit s are plotted in red. The black squares in the left plot are our H2 measurements and 
the empirically determined Lya fluxes presented by Schindhclm et al. (2012b). The green squares are an extrapolation of the -F(H2) vs. 
-F(Lya) relationship. 



to observe a decrease in the H2 content as a function 
of system age. The evolution of the H 2 luminosity is 
plotted in the top panel of Figure 7. The H 2 detections 
are shown in black and the non-detections in red. The 
total H2 luminosity decreases wit h time , in ag reement 
with the findings of llnglebv et all (|2009l I2011al) . How- 
ever, the large scatter in the relation prevents one from 
inferring a characteristic timescale for this decrease. The 
dominant sources of uncertainty in this relation are the 
correlated errors on the age and extinction. The ob- 
served decrease in L(H2) as a function of time does not 
necessarily correspond to a monotonic decrease in the 
H2 content. We expect that the Lya luminosity should 
track the accretion rate (e.g. Fang et al. 2009 demon- 
strate the correlation between accretion and H B aimer 
emissi on), therefore as the accretion rates decline ove r 
time (jArmitage et alj|2003t rSicilia-Aguilar et all 120101 ), 
the associated reduction in Lya emission could mimic 
the effect of H2 dissipation with age. The true situation 
is likely a combination of these effects: gas disk dissipa- 
tion leads to lower accretion rates that in turn produce 
fewer Lya photons to pump the observed fluorescence. 

The lower panel in Figure 7 shows the ratio of the H2 
flux and the Lya flux as a function of system age for the 
14 ob j ects with computed Lya fluxes (jSchindhelm et afl 
l2012bl) . This plot shows the fraction of the stellar + 
shock Lya that is reprocessed by H2. The actual repro- 
cessing factor is i qF(H.2)/F(Lya) > where 77 is a correction 
factor to account for anisotropies in the system geome- 



try and fluorescent radia tive transfer (|Wood et al.ll2002l : 
IWood fc Karov ska 200l . In the case of isotropic emis- 
sion, 77 is simply the geometric filling fraction of the H2, 
as seen by the Lya photons. For this comparison (and 
for the computation of L(Lya) in Figure 6), we have as- 
sumed r\ = 1. For sources where the majority of the H2 
resides in a flattened disk, r\ is most likely less than one 
(Herczeg et al. 2004 found r\ w 0.25 in model fits to the 
spectrum of TW Hya), while for sou rces with a s ignifi - 
cant outflow component r\ may be - 1. lYang et all (pOlTl) 
and iFrance et al.l (|2012| ) present analyses of H2 absorp- 
tion lines imposed on the Lya profiles of the CTTSs 
V4046 Sgr, DF Tau, and AA Tau, systems with both 
high and low inclinations. This implies that at least 
some portion of the H2 in these systems has a Lya cov- 
ering fraction of near unity, suggesting that 77 ~ 1 even 
in some disk-dominated systems. Furthermore, in cases 
where Lya has been scattered out of our line of sight 
or self-absorption redistributes the fluorescence to the 
highe r v levels we use to determine FQI2), r\ can be 
> 1 (IWood et all 120021 ). In the absence of' a more so- 
phisticated radiative transfer treatment of each system 
individually, we assume r\ — 1 for the present analysis. 

The total Lya flux is the full, unabsorbed Lya pro- 
file as it is emitted from the immediate stellar environ- 
ment. Figure 7 shows that the ratio of the total emit- 
ted Lya flux that is reprocessed by H2 is 6.2 ± 2.1%. 
A more meaningful measure of the degree of H2 re- 
processing is the ratio of incident Lya that arrives at 
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Molecular and Atomic Emission Line Parameters. 
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a Average FWHM and average H2 radius (see §3.2) were calculated from four lines of the H2 B 1 ^^ 
H2 luminosity is the sum of the 12 progressions measured in this work (Table 2). 
b Inner H 2 radius, denned a s . R in = GM , (sin(i )/(1. 7 X HWHM[ U |)) 2 . 
c CO line-widths taken from I Sal vk et al.l l!2011aD and lBast et al.l j 201 IB . 

H 1 Lya luminosities with error bars were calculated from the Lya fluxes presented by [Schindhclm et al. (2012b]). Values without error bars 1 
extrapolated from the relationship between -F(H2) and F(Lya). 
e Measurement error on the C IV flux is taken as 5%. 

^ Targets where molecular outflows may contribute to the observed H2 line-width. 

s Strong molecular outflows in RW Aur contaminate the Gaussian fitting (sec Figures 3b and 8). 

^ Targets below the double line do not show measurable H 2 emission in their far-UV spectra. 



the molecular material. The Lya profile will experience 
some degree of absorption in the cir cumstellar environ- 
ment prior to reaching the molecules ([Wood fc Karovskal 
I2004t iHerczeg et al.ll2004t iSchindhelm et al.ll2012al) . Us- 
ing the incident Lya profile observed by the H2, we 
find that the reprocessing fraction (^(H^)/ ^F(Lya)) is 
11.5 ± 1.8%. Therefore, modulo the factor of rj, we 
infer that H2 is capable of reprocessing > 10% of the 
incident Lya flux. This is interesting because the H 2 
a) will isotropically redistribute the Lya photons and 
6) represents a means for transferring Lya photons out 
of the Lya line-core and redistributing them across the 
!000 ^ A < 1650 A bandpass. Additionally, H 2 scat- 
tering is a means for redirecting Lya photons initially 
on a grazing incidence trajectory, increasing the far-UV 
radiation penetration depth. This will significantly al- 
ter the radiative transfer of this fraction of the Lya en- 
ergy as it diffuses outward and downward thro ugh the 
disk dFogel et al.l f2CTlll IBethell fc Bergh] l27)Tl . The 
transfer of these H2-redistributed Lya photons will be 



regulated by circumstellar grains, and this process will 
add power to discrete wavelengths in the far-UV spec- 
trum that propagates towards the disk midplane, possi- 
bly perturbing disk chemistry in regions of active planet 
formation. 

4.2. Spatial Distribution of H2 
4.2.1. H 2 Outflows 

A single Gaussian emission line describes many of the 
observed velocity profiles, however several targets show 
evidence for molecular outflows i n the form o f addi tional 
rcd/blue-shifted H2 emission. iBeck et al.l (2008) pre- 
sented an outflow-selected sample of CTTSs that dis- 
play spatially extended near-IR royibra tional H2 spectra. 
Furthermore, fPontoppi dan et aLI f|201 lh have found slow 
(5 - 10 km s _1 ), weakly collimated molecular winds to 
be common in CO spectra of CTTSs, and varying contri- 
butions from these winds/outflows could be responsible 
for the blue-shifted H2 that is observed towards some 
systems. The COS observations of DF Tau, DK Tau, 
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Fig. 7. — We plot the time evolution of the H2 luminosity 
from the protoplanetary environment (top). The H2 detections 
are shown in black and the non-detections in red. The lower plot 
presents a measure of the H2 reprocessing of the stellar Lyo emis- 
sion, expressed as the H2 flux, d ivided by the empirically deter- 
mined photo-exciting Lya fluxes (Schindhclm ct al. 2012b). r] is a 
correction for radiative transfer effects and the anisotropy of the 
system, set here to r] = 1 (see §4.1). The average H2 reprocessing 
fraction of the total Lya profile is 6.2 ± 2.1 %, assuming 7] = 1, 
shown with solid (mean) and dashed (1 ct) gray lines in the lower 
panel. 



HN Tau, LkCal5, IP Tau, RECX-15, RU Lupi, and RW 
Aur all display H2 line wings extending to the blue of 
the stellar radial velocity. Interestingly, half of these ob- 
jects are known binaries (DF Tau, DK Tau, HN tau, and 
RW Aur), and this may indicate that interactions with 
a companion star contribute to the generation of molec- 
ular outflows in CTTSs. Clearly, more data is required 
to test this connection. 

Of the targets displaying H 2 emission with extended 
blue wings, DF Tau, DK Tau, HN Tau, IP Tau, RECX- 
15, and RW Aur also have reasonably high S/N [0,1] 
progression emission lines in their spectra. The [0,1] line- 
widths are statistically narrower than those from [1,7] 
and [1,4], also suggesting that outflow contributions to 
the [1,4] and [1,7] progressions are present in these tar- 
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Fig. 8. — COS H2 line profiles of RW Aur display a complicated 
kinematic behavior. Stellar variability compromises radial velocity 
measurements for R W Aur. We display +14 km s — 1 , the canonical 
stellar radial velocity (Hartmann ct al. 1986), as the dashed orange 
line. The H2 line profile i s dominated by the redshifted outflow 
lobe (Mclnikov ct al. 2009) with an outflow velocity of m 100 km 



gets (|Walter et all 120031) . RW Aur shows the strongest 
outflow emission in the survey, with both red- and blue- 
shifted emission observed. It is not clear that there 
is a narrow disk H2 component present in this source. 
We show a blow-up of bright lines from the [1,7] and 
[1,4] progressions in RW Aur in Figure 8. RW Aur is 
known to have a bipolar outflow, with the red compo- 
nent (v out ~ +100 k m s" 1 ) being brighter and higher 
densi ty than the blue (jHirth et al.lll9 94; Mcl nikov et al.l 
2009). Due to the COS aperture vignetting function, 
the observed spectra will be dominated by the inner 
w 1" of the RW Aur jet. Figure 8 shows that the 
H2 emission peaks at an observed velocity 80 - 110 
km s _1 to the red of the stellar velocity (+14 km s _1 ; 
Hartmann et al. 1986), suggesting that the molecular 
emission arises in material that is approximately cospa- 
tial with the forbidden atomic line (e.g., [S II] A6731 
A) emission (IWoitas et all 120021: IMelnikov etaLl 120091: 
lHartigan & Hillenbrand! I2009D . The near-IR H 2 outflow 
from RW Aur is centered near ~ +44 km s _1 (jBeck et al.l 
2008), significantly bluer than the peak of the far-UV H 2 
velocity profile. It is not clear if this indicates a differ- 
ence in the physical structure of the UV and IR-emitting 
H 2 , or can be attributed to blending of low- and high- 
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Fig. 9. — A comparison of archival STIS spectra of H2 emis- 
sion lines from RU Lupi (green) and new observations from COS 
(black). The stellar radial velocity is indic ated with the orange 
dashed line and the outflow velocities from Hcrczcg ct al. (2006) 
are shown as the dash-dotted blue lines. Strong lines from the [1,4] 
and [1,7] progressions are shown in the top two panels, and the 
bottom panel shows the two-component spectral fit for the STIS 
observations in blue {Hcrczcg ct al. 2006). This profile, convolved 
with the 1490 A COS LSF is shown as the red dash-dotted line 
and the observed COS data are shown in gray. The COS data are 
consistent with the stellar radial velocity and do not display asym- 
metric outflow profiles as strong as was observed by STIS. This 
difference may be due to a physical change in the inner regions 
of RU Lupi or could be explained by angular extension along the 
dispersion axis in the COS data (§4.2.1). 



velocity gas in the lower spectral resolution near-IR data. 
For the remaining targets, the blue- wings are relatively 
weak, typically only perturbations from the narrower, 
presumably disk-dominated H2 velocity profiles. Further 
study using coadded spectra from several progressions 
would be useful for clarifying the outflow contribution 
and structure in these targets. 

The observed H 2 velocity profiles of RU Lupi are puz- 
zling. Using higher spectral resolu tion observations from 
STIS. lHerczeg et all (|2005U I2006D found that essentially 
all of the fluorescent H2 was contained in two blue-shifted 
components, - 12 and - 30 km s _1 relative to the ra- 
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Fig. 10. — The time evolution of the average radial position of the 
H2 emission in our sample, as defined by the Gaussian line- widths 
of the [1,7] (pumped through the (1 - 2) R(6) absorption line) and 
[1,4] (pumped through the (1 - 2) P(5) absorption line) progres- 
sions (see §3.2 for a discussion). Only targets with disk inclinations 
i > 15° were considered. The red dashed line represents the em- 
pirical relation between the average radius of the H2 emitting disk 
and the system age (Equation 3) for the [1,7] and [1,4] progressions. 
The two H2 progressions show a similar behavior with system age. 
The Spearman rank correlation coefficient for increasing molecular 
radius with system age for both H2 progressions is f» 0.52. The blue 
diamonds indi cate targets with inclinations derived from sub-m m 
dust continua (Andrews & Williams 2007; Andrews ct al. 2011). 



dial velocity of the star. While a blue-wing is apparent 
in the COS observations, the emission is well-fit by a 
single component at the radial velocity of the star. In 
Figure 9, we compare the COS observations (obtained in 
2011 using the 2.5" diameter Primary Science Aperture) 
with the STIS observations (obtained in 2000 using the 
0.2" x 0.06" slit). The apparent line-center velocity shift 
is « 15 - 20 km s _1 , larger than the zero-point calibra- 
tion of the COS wavelength solution^. Furthermore, the 
line-shape is fundamentally different in a way that can- 
not be explained by resolution differences between the 
two instruments. In the bottom panel of Figure 9, we re- 
construct a model two- Gaus sian profile (solid blu e line) 
using the fit parameters from lHerczeg et ail (|2006f) . That 

14 Calibration of the COS wavelength solution is limited by sys- 
tematic uncertainty in the far-UV detector geometric correction. 
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profil e, convolved with the COS LSF at 1490 A (Kriss 
1201 If) , is displayed as the red dash-dotted curve. CTTSs 
(and in particular RU Lupi; Gahm et al. 2008) are known 
to be time- variable, therefore it is not surprising that the 
peak flux has changed, but the line center and emission 
li ne shapes are not con sistent . 

iHerczeg et al.l (|2006t ) found that the H2 fluorescence in 
RU Lupi was spatially extended, although the small size 
of the STIS aperture makes it difficult to predict the ef- 
fects of spatial extension in the larger COS aperture. An 
aperture offset of > 0.4" would be required for significant 
reduction in the instrumental resolving power, but angu- 
lar extension of the H2 emitting region may be able to 
alter both the velocity centroid and the line- width. In the 
dispersion direction, if we attribute the entire 20 km s _1 
offset to extended emission, this leads to a 0.2" displace- 
ment (at the 24.3 milliarcseconds pixel -1 dispersion- 
direction plate-scale of the G160M mode). We cannot 
constrain the angular extent of the RU Lupi H2 lines 
in the dispersion direction, but we note that the optical 
forbidden line emission ([O I] and [S III) are extended 
to ~ 0.2" at a position angle of ~ 225° (Taka mi et al.l 
2001). If the H 2 emission lines observed in the COS 
spectra are cospatial with the forbidden line emission, 
then this could produce the ~ + 20 km s _1 velocity off- 
set. Therefore, it is possible that spatial extension of 
the H2 emitting gas along this axis could explain the 
differences in line shape between the STIS and COS ob- 
servations. We compared the cross-dispersion profile of 
the two-dimensional spectrogram of the RU Lupi ob- 
servations (over the range AA ~ 1420 - 1450 A) with 
that of the DA white dwarf WD0320-539. The two spec- 
tra were centered to within 0.5 pixels (« 0.05") in the 
cross-dispersion direction and had nearly identical pro- 
file FWHMs (4.5 pixels for WD0320, 4.8 pixels for RU 
Lupi). 

We conclude therefore that either RU Lupi has spatial 
extent along the dispersion axis, or the spectral profile 
differences between the COS and STIS epochs are caused 
by a physical change in the system. RU Lupi may con- 
tinue to be outflow dominated, however the COS obser- 
vations raise the possibility that that RU Lupi was ob- 
served during an episode of strong outflow in 2000 or that 
the geometry has evolved such that disk-illumination by 
Lya contributed more strongly during the 2011 observa- 
tions. Continued spectral monitoring of these loopy line 
profiles would be interesting. 

4.2.2. H 2 Disks 

The high S/N H 2 lines in the majority of our sample 
targets can be adequately described by a single Gaussian 
emission component at the stellar radial velocity. Bear- 
ing in mind the added uncertainty introduced by blue- 
shifted H2 in some targets, we conclude that most of the 
observed fluorescent emission originates in a disk. This 
is in c ontrast with the near-IR H2 study of Bec k et all 
(2008); however, that study targeted stars known to have 
strong outflows, which is not the case with the sample 
presented here. Using the velocity-resolved H2 line pro- 
files, (B.H2) was calculated for the [1,7] and [1,4] progres- 
sions. The radial distributions as a function of system 
age are displayed in Figure 10. The general result is that 
the average location of the emitting H2 gas in the inner 
disk moves outward from ^0.25 AU to ~ 2 AU as the 
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Fig. 11. — A com parison of the emis s ion line FWHMs of H2 
(this work) and CO ISalvk et al.ll2011al: Bast et all 12011). The 
squares are the Gaussian line- widths of the [1,7] progression and 
the stars are the Gaussian line- widths of the [1,4] progression. The 
dashed line represents the expected correlation if the H2 and CO 
line- widths are the same. 
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disk traced by H 2 (this work) and CO ISalvk et al.ll2011al ). The 
dashed line represents the expected correlation if the H2 and CO 
emission are co-spatial. Ri n is defined as the Keplerian radius cor- 
responding to an orbital velocity of 1.7 X HWHM of the Gaussian 
line fits (§3.2). Rin{^2) is computed from the [1,7] progression 
for this comparison. Except in the case of SU Aur, i?i n (H2) and 
-Ri n (CO) are the same to within a factor of three, suggesting that 
the observed CO and H2 orbit at similar radii. 



system evolves from 1-10 Myr. 

The relationship between the system age and the av- 
erage H2 radius can be characterized by the simple em- 
pirical formula 

log 10 (R(H 2 )) m = 0.81 log 10 (Age) - 5.36 (3) 
where (Rh2)[i,7] is in AU and the age is in years. The 
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coefficients in Equation 3 [-5.36 ± 1.91, 0.81 ± 0.30] 
are computed from a y 2 minimization of a linear func- 
tion of logio(Age) and logio(i?//2)[i,7]- The Spearman 
rank correlation coefficient for increasing molecular ra- 
dius with system age is 0.52 with a deviation from zero 
of 1.1 x 10~ 2 , meaning that there is a relatively low prob- 
ability that logio(Age) and logio(-Rf/2)[i.7] are uncorre- 
cted. Targets with disk inclinati ons derived from sub- 
mm dust continuum obser vations ([Andrews fc Williams! 
120071: lAndrews et al J 1201 lh are plotted as the blue dia- 
monds in Figure 10. There is a large spread in the distri- 
bution of {Rh2) with age. We suggest that in addition 
to contamination by outflows, uncertainties in the stel- 
lar mass, inclination angle, and stellar age all contribute 
to the dispersion in this correlation. With the exception 
of CS Cha (§4.3), the H2 in all disks is concentrated at 
a < 3 AU. 

The weak trend towards larger radii as a function of age 
suggests a scenario where the average molecular emission 
radius moves to beyond pa 1 AU in about 10 Myr. It is 
not immediately clear how to interpret this result, and we 
remind the reader of the ingredients necessary to produce 
H2 fluorescence in a disk. The H2 opacity must be high 
enough to absorb a significant number of Lya photons, 
requiring both appreciable column densities and a suf- 
ficient population of H2 in excited ro-vibrational states. 
The latter requires a hot (T(H 2 ) > 2000 K) molecular 
layer, possibly with a contribution by intense illumina- 
tion from the A < 112 A continuum from the central star 
+ ac cretion shocks (jNomura et al.l I2007t iFrance et al.l 
2012). The second major requirement for observable flu- 
orescence is a geometry where the disk subtends a sub- 
stantial angular cross-section of the Lya emitting area, 
that is, the disk must be sufficiently flared in order to in- 
tercept enough Lya photons to excite the observed emis- 
sion. The observation of larger H2 radii as a function 
of time suggests that one or more of these criteria are 
not being met in the inner pa 1 AU in the more evolved 
systems. It may be that the hot H2 is in the process of 
dissipating from this region, possibly due to dynamical 
clearing by a protoplanet, enhanced H2 dissociation as 
this region is less shielded by grains, or photoevapora- 
tion by energetic radiation from the central star. Alter- 
natively, this result may indicate that the flaring angle 
in the inner disk is decreasing across our sample, sug- 
gesting an evolution of the vertical structure of the disk 
on timescales of a few Myr. Improved stellar masses and 
ages, and larger samples of well-determined disk inclina- 
tions would allow better characterization of the relation 
(or lack thereof) between molecular radius and age, en- 
abling a better understanding of the evolution of inner 
gas disks. 

4.2.3. Comparison with near-IR H2 Emission 

Quadrupole rovibrational line emission from H2 (most 
notably the (1 - 0) S(l) A2. 1 /xm line) has been detected 
around several CTTSs feary et a l. 2003; I toh et al.ll200l 
ICarmona et~aT1l2007t IBarv et al.ll2008D . While the num- 
ber of objects available for direct comparison is small, 
it seems that the Lya-pumped H2 emission arises in- 
terior to the near-IR emission lines (although see §4.3 
for the case of the transitional disk CS Cha). Typical 
emission line widths for the near-IR H2 sample presented 



in IBarv et all (12008! ) are FWHM < 20 km s" 1 , and are 
thought to originate at radial dis t ances of a few to a few 
tens of AU. Similarly. IBarv et al.l (|2003ft observed (1 - 0) 
S(l) emission in LkCal5 with FWHM < 14 km s^ 1 and 
suggested that the emitting gas was located between 10 - 
30 AU from the star. This is a factor of ~ 4 narrower 
than the FWHMm 7 ] — 53 ± 3 km s _1 that we observe 
in the LkCal5 UV spectrum. 

RECX-15 is the only disk in the 77 Cha region to 
emit a measurable flux of Hg in the (1 - 0) S(l) 
line ()Ramsav Howat fc Greaves! [2007h . The near-IR H 2 
line widths (18 ± 1.2 km s _1 ) are a factor of pa 2.3 
smaller than measured in the COS spectra of RECX- 
15 (41 ± 4 km s _1 ), although we note that outflows do 
contribute to the RECX-15 ultraviolet H2 spectra. The 
near-IR H2 emission from disks is typically interpreted as 
disk surface gas excited by energetic radiation. Includ- 
ing t he effects of gr a in gro wn, UV, and X-ray illumina- 
tion, iNomura et al.l (|2007l) have demonstrated that gas 
temperatures in the range 1500 - 3000 K can be main- 
tained in the disk surface to radial distances of ~ 10 
AU from the central star, therefore the excitation con- 
ditions necessary to both produce near-IR emission and 
to enable Lya-pumping appear to exist to at least this 
radius. However, the fluorescent ultraviolet emission is 
dominated by H2 nearest to the source of Lya photons. 

4.2.4. Comparison with CO Emission 

Comparing the line-widths of different molecular 
species can provide an observational constraint on the 
composition and physical structure of inner gas disks. 
Spectral observations of the 4.7 fim fundamental band 
CO emission are a widely u sed tracer for this mate- 
rial (jSalvk et al.l 120081 120091 ) . and understanding the 
molecular structure and the degree to which various spec- 
tral diagnostics trace the same gas are useful towards 
a more complete picture of the planet-forming regions 
around CTTSs. Figure 11 shows a comparison of the 
Gaussian FWHMs of CO and H2 for the subsample of our 
targets that have been observe d by high-resolutio n mid- 
IR spectrographs (|Salvk et alJl2011allBast et al.ll2011[ ). 
The emission line- widths from the [1,7] and [1,4] H2 
progressions are self-consistent in all cases and approx- 
imately equal the CO line-widths up to FWHM pa 60 
km s -1 . At larger CO line- widths, the H 2 FWHMs do 
not exceed 70 - 80 km s _1 , which may indicate a phys- 
ical boundary condition inside of which H2 is subject to 
collisional and/or photodissociation. 

The inner radii of the H2 and CO disks can also be di- 
rectly compared. The calculated H2 inner radii are pre- 
sented in Table 3, and Figure 12 compares th e ultraviolet 
H 2 (U V-H 2 ) and infrared CO (IR-CO) radii (|Salvk et al.1 
l2011al ): the dashed horizontal line represents a one-to-one 
relation. With one exception, R in (R 2 ) and i?i„(CO) are 
the same to within a factor of three. The agreement be- 
tween the two molecules is rather remarkable given that 
we are comparing different species, excited by different 
mechanisms (photo-excitation vs. collisional excitation), 
observed at different epochs in different wavebands. The 
notable exception is SU Aur, whose i?m(H2) is approxi- 
mately an order magnitude larger than its corresponding 
-Rirt(CO), possibly due to H2 emission from nebulosity 
associated with this star. The agreement between the 
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IR-CO emission and the UV-H2 emission is also interest- 
ing in light of the recent discovery of large amounts of CO 
emission in th e UV spectra of CTTSs ( UV-CO; France 
et al. 2011b). iSchindhelm et al.1 (|2012al ) present an ini- 
tial survey of this emission, showing that the line-widths 
of UV-CO are systematically narrower than those of 
UV-H2. Their interpretation favors a picture where the 
UV-CO originates in a cooler molecular layer (T rot (CO) 
~ 500 K) at larger semi-major axes (a > 2 AU) than 
both the UV-H2 and the IR-CO, consistent with the re- 
sults presented here. 

4.3. H2 in Transitional Disks 

The majority of H2 emission in the targets in our sam- 
ple originates from a < 3 AU. A notable outlier from 
the average H2 radii presented in Table 3 is CS Cha, 
which lies at significantly larger (Rh2)[i,7] than the rest 
of the disks studied here. The [1,7] line- width of CS 
Cha is 18 ± 7 km s _1 , which makes it the only unre- 
solved moderate inclination targeiP^I. and allows us to 
place a lower limit on the average H2 emission radius, 
(■Rfl"2)[i,7] > 9 AU. CS Cha is a transitional disk, showing 
the largest mid-IR spectral slope in our sample, ni3_3i 
= 2.89 (Furlan et al. 2009; see Table 3). Modeling of 
the Spitzer-IRS mid-IR spectrum of CS Cha reveals a 
trun cation of the inner dis k dust distribution at a « 43 
AU ([Espaillat et al.ll2007bT) . possibly the result of a dy- 
namical interaction with a companion star (Guenther 
et al. 2007, but see also Espaillat et al. 2011). The 
new HST data presented h ere confirms the presence of 
molecular gas in the system ([Barv et al.ll2008l) , and while 
we cannot rule out the possibility that this gas is co- 
spatial with the dust at a > 40 AU, this would imply 
that hot molecular gas (T(H 2 ) > 2000 K) exists at large 
radial distances from the star. Therefore, because only 
(A(H2) > 10 18 cm -2 ) is required to produce detectable 
fluorescence, this emission may originate in the tenuous 
molecular material in the disk gap, or if Lya photons can 
propagate through the gap, this emission may arise from 
the edge of the directly exposed wall at ~ 43 AU. 

The H2 emitting material we observe in CS Cha may be 
ph ysically associat e d with the [Ne II] emission observed 
by lEspaillat et al.l (|2007b[ ) , but that cannot be conclu- 
sively determined from the available observations. We 
can rule out an origin in the optically thin dust disk in- 
side of 1 AU, as our limit on the H2 inner disk radius in 
CS Cha is R in (R 2 ) > 3 AU. 

CS Cha is one example of a generic property of our 
sample: H2 is common in the inner regions of accreting 
transitional disks (ni 3 _ 3 i > 0.5). In addition to CS Cha, 
our sample includes the well-studied transitional systems 
DM Tau, CM Aur, UX Tau A, LkCal5, HP 135344B 
and TW Hya ([Calvet et al.l [20051 l20l)l lEspaillat et aD 
[20073). H 2 emission is found to originate inside the 
dust hole in these systems, (Rh2) < Rdust, consistent 
with the origin of the IR-CO gas in t he inne r regions of 
transitional disks ([Salvk et al.l 120091 l2011aft and previ- 
ous observations of n ear-IR rovibrational emission from 
H 2 (IBarv et al.ll2"008h . 



5. SUMMARY 

We have presented the most sensitive survey of H2 in 
protoplanetary environments to date. The majority of 
this work was made possible by the combination of large 
effective area and low instrumental background at mod- 
erate spectral resolution provided by the i?ST-Cosmic 
Origins Spectrograph. We have used this survey to mea- 
sure the time evolution of both the spatial distribution 
and the luminosity of H2 in young, low-mass disks for 
the first time. Below we summarize the primary results 
of this work: 

1. We obtained far-UV spectra of 34 T Tauri stars: 
27 accreting CTTSs and 7 non-accreting WTTSs. 
Of these, 100% of the accreting sources display a 
measureable amount of H2 emission, providing di- 
rect evidence for the interaction of a strong Lya 
radiation field with the molecular disk surface. 

2. We found that the H2 luminosity is well correlated 
with the Lya and C IV luminosity, consistent with 
a scenario where gas-rich disks fuel larger accretion 
rates that produce energetic radiation. 

3. The H2 luminosity is observed to decline with age, 
although H2-rich systems persist to ages ~ 10 Myr. 

4. We measured resolved H 2 line profiles of 23 tar- 
gets with inclination angles > 15° and found that 
these line profiles are reasonably well fit by a single 
Gaussian component at or near the stellar radial 
velocity. Assuming a disk origin for these targets, 
we used the line- widths to constrain the spatial dis- 
tribution of the emitting molecules to a < 3 AU in 
most cases. 

5. The inner radii of H 2 disks are roughly consistent 
with those of CO disks (measured from A ~ 5 fim 
spectroscopy). 

6. Strong H2 emission is observed at a > 0.2 AU in a 
subsample of transitional disks (n.13-31 > 0.5). 
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